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We consider Double Quintessence models for which the Dark Energy sector consists of two cou- 
pled scalar fields. We study in particular the possibility to have a transient acceleration in these 
models. In both Double Quintessence models studied here, it is shown that if acceleration occurs, 
it is necessarily transient. We consider also the possibility to have transient acceleration in two 
one-field models, the Albrecht-Skordis model and the pure exponential. Using separate conservative 
constraints (marginalizing over the other parameters) on the effective equation of state We}}, the 
relative density of the Dark Energy f2Q,o and the present age of the universe, we construct scenarios 
with a transient acceleration that has already ended at the present time, and even with no accel- 
eration at all, but a less conservative analysis using the CMB data rules out the last possibility. 
The scenario with a transient acceleration ended by today, can be implemented for the range of 
cosmological parameters Slm.o ^ 0.35 and h < 0.68. 

PACS numbers: 04.62.-|-v, 98.80.Cq 

I. INTRODUCTION 

The release of type la Supernovae data independently by two groups, the Supernovae Cosmology Project and 
the High-z Survey Project jlj, confirmed in more recent work 0|, indicating that our Universe might be presently 
accelerating, has profound implications on the current paradigm in cosmology. It was found that, assuming flatness, 
the best- fit Universe with a cosmological constant A is given by the set of cosmological relative densities (riA,Oj f^m.o) = 
(0.72,0.28). Here, m stands for usual pressureless matter including (cold) dark matter and baryonic matter. These 
data are also in surprising agreement with the location of the first acoustic (Doppler) peak of the Cosmological 
Microwave Background temperature anisotropy multipoles. If these data are confirmed in the future, they imply 
a radical departure from usual textbooks Friedman-Lemaitre- Robertson- Walker (FLRW) cosmology Indeed, a 
perfect isotropic fluid cannot lead to accelerated expansion unless it has a sufficiently negative pressure while the data 
seem to imply that such a kind of smoothly distributed matter, called Dark Energy, constitute about two thirds of 
the whole matter budget of our Universe. 

Obviously, a pure comological constant A could be responsible for this acceleration. However, its amplitude has to 
be exceedingly small, about 123 orders of magnitude too small in order to be explained in a "natural" way. While 
this possibility is actually in good agreement with observations and attractive as it seems to make many pieces of our 
present understanding of structure formation fit into a consistent picture, in view of the above mentioned theoretical 
problems research has focused on other, more elaborate. Dark Energy candidates that would mimic a pure cosmological 
constant. 

A A-term is equivalent to a perfect isotropic fluid with constant equation of state pA = — PA, equivalently Wf^ = —1. 
Provided enough of the total energy density is stored in this component, the expansion will be accelerated. A first 
generalization could be some phenomenological scaling Dark Energy with constant equation of state and w < 
Some interesting insight can be gained by considering most of the universe energy stored in such a component (see 
e.g. i). 

Another more elaborate alternative to a pure cosmological constant A is some effective, slowling varying, cosmo- 
logical constant term that will start driving the universe expansion at low redshifts. The prominent candidate in this 
respect is some minimally coupled scalar field often termed quintessence, slowly rolling down its potential so that 
it has a negative pressure, p$ = it;$p$, w,^ ~5 [SIHQ, 1^. Later on tracking solutions were introduced 9], a sub- 
stantial improvement with respect to the initial conditions problem, however without solving the cosmic coincidence 
problem. Actually, this is precisely the mechanism which drives the inflationary stage in the Early Universe and here 
too, one can consider several potentials and investigate how well they fit the observational data. If a scenario of 
this kind is the correct one, then it is possible in principle to reconstruct (the relevant part of) its potential V{'^) 
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and the corresponding equation of state characterized by w$ , using luminosity distance measurements in function of 
redshift z. This procedure that can be extended to more elaborate models of gravity like so-called generalized, or 
extended, quintessence models in the framework of scalar- tensor theories of gravity llflllll . Note that the latter case 
can represent a possible physical realization of so-called phantom energy, with w < — 1 lCr| and it has been intensively 
investigated (see e.g. One can also consider other modifications to the theory of gravity or to the coupling of 

scalar fields (see e.g. 

A Chaplygin gas is another example of more exotic Dark Energy candidates that was proposed, in the course of the 
universe expansion, this gas undergoes a transition from a dust-like to a cosmological constant-like equation of state 

Actually, as the Dark Energy sector is still unknown, it is interesting to explore all possible models that can pass 
successfully the observational constraints and we want to investigate Double Quintessence models for which this sector 
consists of two coupled scalar fields. Two scalar fields models give rise to exchange of kinetic energy between the 
two fields and produce a non-trivial time evolution of the Dark Er iergy equa tion of state. Some models of two fields 
quintessence have been introduced and studied in the literature We show with simple two-fields models 

how the introduction of an auxiliary field can bring accelerated expansion of our Universe to an end. Though our 
models have natural parameters, the cosmological coincidence however is not solved here. We note that many-fields 
infiationary models can have signatures |19|| which distinguish them from single-field inflationary models, in the first 
place the possibility to have a characteristic scale in the primordial fluctuations spectrum (see e.g. |20]). In this 
respect, the situation is here more contrasted and we will return to this point in the Conclusion. 

The existence of an event horizon in the case of a de Sitter phase is an obstacle to the implementation of string 
theory because the S-matrix formulation is no longer possible, hence eternal acceleration, leading to a de Sitter space 
in the asymptotic future is problematic plj . Therefore, transient acceleration of our universe is certainly a wellcome 
feature in this respect. Even more, we will show that we can have, in a fairly natural way, a transient acceleration 
which has already ended at the present time and still in accordance with the observations. We explore the window in 
parameter space for which this intriguing scenario is realized. It is even possible, though this is a marginal possibility, 
to construct scenarios where there is no acceleration at all, neither at the present time nor in the past of our Universe. 
Of course, agreement with the observational data still require that the Dark Energy sector dominate from some time 
on and has an effective equation of state satisfying w < 0. We will see that in our models the quantity w, which 
will be called eos (equation of state) parameter, is not constant and actually strongly varying on the redshift interval 
< z < 2 relevant for (past) luminosity distance measurements. Actually, it was already found earlier that luminosity- 
distance measurements are rather insensitive to large variations of the eos parameter at low redshifts 0, ■ This is 
why, very interestingly, in our Double Quintessence models, it is possible to construct scenarios where the standard 
interpretation of a presently accelerating Universe can be challenged. 

In Section II, we give the basic equations of our system and introduce the relevant quantities and notations. In 
Section III, we list the observational constraints on our models and the associated quantities. In Section IV, our 
models are presented and investigated numerically in Section V and regions in parameter space leading to viable 
models are explored. We will take SttG = M'"^ = 1 (where Mp is the reduced Planck mass) and ft = c = 1 in the 
following. 



II. COSMOLOGICAL EVOLUTION 

Let us introduce now the background evolution equations of our system. The equations of motion for a spatially flat 
FLRW cosmology with two coupled homogeneous scalar fields <& and 5", and Hubble parameter i7, are conveniently 
written in the following way: 

Ph = -iH-fbpb (1) 
$ = -iH^-d^V (2) 
* = -iHi!-d^V (3) 
H = -47rG'(7„p,„ + 7rPr + IqPq) , (4) 



subject to the constraint equation: 



SttG k 

H = —^(P'm+ Pr + Pq) - ■ (5) 



Here a dot denotes a derivative with respect to the cosmic time t, the subscript b refers to the dominant background 
quantity, either dust (m) or radiation (r) while Q refers to the Dark Energy sector, here the two quintessence scalar 



3 



fields. We have further 

Pb = [lb ~l)pb , (6) 
with 7m — 1 and 7r = | where, for any component i, we have introduced for convenience the quantity 

-fi = l + Wi . (7) 
FinaUy the quintessence fields with potential V have the following energy density and pressure: 

PQ^l^' + l^^ + V{<S>,'f) (8) 

PQ = ^*' + ^*'-n*,*) (9) 
with pQ = (7(3 — l)pQ. It is convenient to define the following new variables : 



8ttG $ ^ /SvrG * ^ 8ttG r- , , 

The above equations are then written in the following way : 

$' = (87rG)"^\/6 X$ (11) 
^' = (87rG)"^%/6 X* (12) 

= -3X$ - y|(87rG)-^ X2 a^lny + ^ X*^^ (13) 

= -3X* - y|(8^G)-^ X2 W + ^ (14) 

-(87rG)"^Xv(X$9$lny + X*a^.lny) + - XyF (15) 

with 

i^(X*, X*, Xv, iV) = {\-Xl-Xl- Xl) (IjuPhiUiPl] + 2(X| + Xl) (16) 

where a prime denotes a derivative with respect to the quantity N ^ the number of e-folds with respect to the present 
time, 

iV = In— , (17) 
ao 

and we have also H = N . It is straightforward to relate the quantity N to the rcdshift z 

l + z^e'^ . (18) 

The fields are expressed in units of the reduced Planck mass and all quantities above are dimensionless. Moreover we 
have the following standard evolution for matter and radiation densities 

Ptu = Pm,o exp(-3A^) , (19) 

Pr = p^,oexp(-47V) , (20) 

and more generally for constant wx 

px = pxfl exp(-3[l + ivx]N) . (21) 

Hereafter and in p7|) . the subscript refers to the value of any quantity today and the subscript i refers to its value 
at some initial time t, . 
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III. OBSERVATIONAL CONSTRAINTS 



We consider now all the relevant quantities pertaining to our system which will enter the observational contraints. 
The relative energy density for matter, radiation and quintessence, f2,„, and flq, where we have for each 
component labelled by the subscript i, fi^ = ^ with p^, the critical density, are given by 



"m,0 + "r,0 e 



— - Xq, - - Xy)— — j , (23) 

Qq = Xl + Xl+ X2 . (24) 



The equation of state (eos) parameter wq for the Dark Energy (Double Quintessence) sector, and the effective eos 
parameter Weff [l^l read 

_ PQ - Xy 

^0 = — = ^2 , y2 , v2 ' (25) 



PQ Xi+X^+ X 



V 



^ C da'wQ{a')nQ{a') /"^ dN' e^' {Xj + Xl - Xl) 
"^^^ C da'n^ia') dN'e-'iXl + X| + X^) ' ^''^ 



The case where Dark Energy consists of one field is straightforwardly recovered from the above equations. Finally, 
the deceleration parameter q, the Hubble-parameter-free luminosity distance Dl and the age of the Universe to are 
respectively given in terms of A'^ by : 

i 

l{l + Qr + iWQflQ) , (28) 



2 

Dl = HodL ^{l + z) I dz'-^ = e-'' I dN'e-"' , (29) 

Jo J^V- ) JN 

+ 00 TT TT 







with 



Hq J 1 — X| — X^ — Xy 



(31) 



Let us consider two universes with the same "history", in the sense that they share the quantity h{z) = while 
they differ in the present value Hq. As can be seen immediately from ()29|) . the quantity Dl depends only on h{z) and 
will therefore be the same for both universes. It will be interesting to compare luminosity distances o?l(z) for such 
universes, we have in particular when comparing two universes with same h(z) but different Hubble constant Hq and 

dL{z;Ho) = dL{z;HQ,A_) = dL{z;HQ,A_) ee ga dL{z;Ho.A) , (32) 

Ho h 

where h = iqo knT/s/Mpc ' '^'^'^P- = loo km/s/Mpc " ^^^^ further for any two universes (dropping the argument z) 

(/ia) = GA— — , (33) 

ttL.A L)l,A 

which shows the relative variation of dL^z; Ho) when varying Ho with respect to some fixed value i?o,A- We will make 
use of (1^ in Figure 10. 

The traditional sign convention for the "deceleration" parameter q gives a positive q for a decelerating universe. A 
cosmological constant A corresponds to the particular case wa = — 1 and J^a = For flat universes, the relative 



5 



energy densities ili satisfy J^i — + + = 1 at all times. In particular, at late times, the energy density 
of radiation can be neglected and H27I28|) gives 

q ~ ir!„, + hlQ{\ + Zwq) - ^(1 + 3u.Q ^q) . (34) 
We see in particular from (|77|l that a decelerated expansion at late times requires 



WQ > -- . (35) 



In this work, we will consider universes where the quantity q (re)changes sign, from negative to positive, i.e. acceler- 
ating universes resuming a decelerated expansion, in some cases even before the present time. We will call henceforth 
tend the time at which the transient accelerated stage ends 

q{t > tend) > . (36) 

For our system, the eos parameter wq is constrained between — 1 < wq < 1 while the energy density of the Dark 
Energy sector (like the energy density of any component) is bounded according to < fig < 1. The parameter q has 
to be negative at the present time if the universe is accelerating today, viz. qo < 0. 
The nucleosynthesis bound is the most stringent one |23| 

J7Q(iV - -23) < 0.045 at 2(7 , (37) 

while we have at last scattering |2^ 

Qq^N 8) < 0.39 at 2cr . (38) 



We adopt the following range 0.6 < /i < 0.8 and we consider universes satisfying the conservative bounds |2l l25l l26l 

0.2 < n,n,o < 0.45 , (39) 

nr^oh^ = 4.3069 X 10~^ (40) 

as well as 

0.55 < r^Q^o < 0.8 . (41) 

In umin, a constant eos parameter wq was considered. The quantity Weff was introduced in in order to 
account for a varying equation of state of the Dark Energy sector. Note that for constant eos parameter one has 
Weff = Wq and it was shown in psj that constraints on wq can be replaced by constraints on Weff so we consider 
the observational constraint : 

Weff < -0.70 . (42) 

The conservative bound ^ is required by the CMB as well as the SNIa and the LSS data pi I2I I27L I2I. l29 t . 

Combining eqs (|35I41|) , decelerated expansion today requires the necessary condition 

WQ,o > -0.606 . (43) 

We note immediately that conditions H42|l and H43|l are incompatible for a constant equation of state, i.e. Weff — 
WQ = constant. Both conditions (|42I43|) can be met today when the equation of state varies strongly at low redshifts 
as we will see with concrete models studied in this work (see Figures r6l8|l . The lower bound H43|l is always necessary 
but only sufficient when Qqa — 0.55. When fig.o ~ 0.8, the condition sufhcient for decelerated expansion is tighter, 
namely wq_q > — ^ ~ —0.416. Clearly, the lower ilp.n, the easier it is to implement decelerated expansion today. 
Finally we impose that the age of the Universe satisfies [33 

to > 13 Gyrs , (44) 

which translates into a constraint on the quantity Hq, or equivalently on h, as can be seen from H30|) . One has 

= 3.0856 X lO^^/i-^s = 9.7776/i~iGyrs , (45) 
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so that the constraint 1)44(1 can be rewritten as 

Hoto > 1.330 h . (46) 

Note that in the past, the inclusion of a cosmological constant was invoked precisely in order to reconcile an "old" 
universe with a high value for Hq while today it is motivated by completely different observations. The bound 144(1 
is conservative and corresponds to the present age of an Einstein-de Sitter universe with h — 0.65. In our simulations 
we will assume equipartition of the energy 

i$2 _ 1^2 _ y. (47) 

initially at the end of the primordial inflationary stage (z^ ~ 10^^), but the initial time could as well be taken at 
nucleosynthesis or at matter-radiation equality |32l | . 



IV. MODELS 



We will be interested in universes where the accelerated expansion is brought to an end, or even does not take 
place. If the universe is accelerating indefinitely it will exhibit an event horizon such that : 

f°° dt 

Duit') = / 4y (48) 
Jt' a{t) 

is finite. In the case of a constant eos parameter wq and a quintessence domination {Qq ^ Qm,^r), we have 

2 

. oc t^<i+™Q> and if wq < -1/3 the integral gHIl is finite. 

It is well-known that the problem of initial conditions for the quintessence field can be considerably alleviated by 
the so-called tracking behaviour whereas the relative density Jig follows the density f2f, of the dominating background 
component, with some specific evolution of wq. We will consider two one- field models and two Double Quintessence 
models. 

Let us stress first what will be relevant in the framework of either one-field, or two-fields quintessence models with 
transient acceleration that we study here. There are several ways to obtain transient acceleration: 

• the tracking behaviour of the scalar potential is not reached until today, namely fields are frozen in the past with 
Wq — —1, ^Q <^ 1, while once the tracking behaviour is reached, ea. l(^ must be satisfied. In Section IV. A, we 
will see that a pure exponential potential can exhibit such a behaviour as studied in |3^ 133 . 135ll3^ . Moreover, in 
[33 it is shown that a FLRW universe filled with a scalar field only can undergo a stage of transient acceleration. 
As we show in Section IV. D, even the quintessence domination can be avoided because of the varying coupling 
constant of the model studied there. 

• the existence of a minimum or a local flatness in the $ direction of the potential can produce an accelerated 
expansion but the minimum has to diseappear in order to make this acceleration transient. In the model of 
Albrecht & Skordis 38] (section IV. B) for which a feature is introduced in a pure exponential potential, the 
parameters have to be fitted in order for the field to roll over the barrier. In model IV. C, the minimum is 
dynamical and disappears. So in these cases, the actual acceleration takes place outside the scaling regime, 
which was reached early on in the universe evolution, when the field approaches his minimum, see section 
IV.B and IV.C. 

• Still another possibility to obtain transient acceleration is to cancel the scalar potential V : this is the case in 
the hybrid inverse power law potential V{^,'^) = \E'2Af'*+"$^" studied in [13. On the contrary, the inverse 
power law potential cannot produce a transient acceleration because wq ^ — 1, rig ^ 1 in the future 
Other models where V cancels can produce transient acceleration like the oscillatory Dark Energy model with 
a double exponential potential V{^) = (ylexp(iA<l') — S cxp(— ^A^)) [s^ and also the power-law potential 
V (X $2", n= 1,2,... 



If y < 0, even a flat universe can undergo a transient acceleration followed by a big crunch |40|,|4ll, and an ever 
decelerating universe can also be considered in a closed FRLW (fc = -1-1) j43|. In the following we will restrict ourselves 
to positive potentials V > and flat space k = 0. In the transient acceleration picture, it is interesting to consider 
whether the acceleration flnishes in a time comparable to H^^, namely t^nd ^ or tend ^ ^o? nruch larger than H^^, 
i.e. tend 3> in, or if the acceleration does not occur at all. 
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A. Pure exponential potential 



We consider a scalar field $ with an exponential potential 

Vi^) = M^e-^* , (49) 

already widely investigated in the past 0, S 0, and motivated in "il] . For this potential an attractor solution 
exists, either a scaling solution such that : 

^Q = ^ and 7Q = 76 if > 3jb , (50) 

or else a scalar field dominated solution : 

Oq = 1 and 7Q = y if < < 87^ . (51) 

The nucleosynthesis bound i^gil MeV) < 0.045 at the 2a level implies from (|50|l that A > 9 during the scaling 
regime. If we want to take advantage of the attractor property of this potential, it is impossible to have reasonnable 
^Q,o without violating the nucleosynthesis bound, and anyway wq would mimic the background component eos, 
WQ = Wm = thereby preventing a past and/or actual acceleration. 

In [33I I33 . 135I I the authors present a way to circumvent these arguments and revive the pure exponential potential, 
using the attractor property in the future instead of the past of the universe by fine tuning the mass scale ^ MpHo 
(in this Section we will sometimes reput Mp for clarity). Contrary to the claim in [s^l, $^ ^ 1^ is not necessary, only 
Vi ~ MpH^ has to be imposed initially. Apart from A and M, two extra values and $i must be specified. We can 
always take $i = by a redefinition of M which corresponds to a rescaling of the problem. 

The universe experiences first a kination regime such that V <^ (x a"^ with wq ~ 1 followed by a regime 
during which V ^ and while Mpm\ ^ V ^ MpH"^, the field gets frozen with wq ~ —1 until now where 
$2 ^ Mpm| ~ 1/ ~ M^H^ allowing the attractor regime to be reached in the near future and we are left with a fixed 
ratio, either fig = 1 and wq = A^/3 — 1 if A^ < 3 or else 57q = 8/A^ and wq = if A^ > 3, in the presence of dust 

miiiiii. 

In contrast to the statement in [s^l, the value A > VS gives a viable model provided the scaling regime is not yet 
reached and it produces intriguing possibilities for the evolution of the universe, namely an acceleration which ends 
before today and even no acceleration at all. 

An acceptable quintessence density ilg > ilg.o ^ 0.55 implies, using (|50|l . A < 2.3. If A < V2, acceleration is 
eternal because the attractor solution (in the future) is characterized by VIq = 1 and wq < —1/3 as seen from (|51|l . 
Therefore transient acceleration implies ^/2 < A < 2.3 and an eternal quintessence density domination (fig > flm). 
Numerical results of section V.A can tighten even more the upper bound. This model has one fine tuned parameter 
M, which has to satisfy AP MpHo. 

It is possible to produce analogous scenarios containing multiple scalar fields with an exponential potential ^5l| of 
the type 

n 

1=1 

In the context of assisted inflation, though each field is unable to support separately an inflationnary stage, together 
they are able to do so. 

But we are interested in a transient acceleration whether already finished or not, hence in non inflationnary solutions 
at late times, which requires the condition 



n 

i=i ' 

because both attractor solutions H50I51|) remain valid with the change 

(54) 
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If each $i = initially, then M has to satisfy M ^ ^JM^Hq with some level of fine tuning. 
Models containing coupled scalar fields with an exponential potential of the type 

= Af^e-SJLiA,*, ^ (55) 
can produce the same scenarios with replaced by X^jli i'^ expressions H50I51|I . So if 

transient acceleration will occur even though for each slope separately Aj < a/2- Again if each $j = initially, then 
M has to satisfy M ~ y/MpHo. 

It is possible to generalize the two last cases '^i^h the potential 

n 

V = M^Y1 e~^^"=^ • (57) 

1=1 

It can be shown that in the presence of a barotropic fluid, the two late-time attractor solutions are, either a scaling 
solution : 



fig = ^ and 7q = if A^ > 37f, , (58) 

or else a scalar field dominated solution 



with 



\2 

1]q = 1 and 7Q = y if < A^ < 87^ . (59) 
1 " 1 

T2 = H v"' • (60) 

For example if each A^^ = A and rrii — m then A,. = \/^A. Thus assisted inflation tends to lower A,, and the coupled 
part to increase it, but when 

2 < ^ 1 , (61) 

the late-time solution produces a decelerated expansion with the fine tuning M ^ ^MpHo if initially each $y = 0. 



B. Albrecht & Skordis potential 

Albrecht and Skordis js^ proposed an interesting model of quintessence. As noted by this model contains 
solutions for which there is a transient acceleration of our universe. We would like to emphasize even more that 
acceleration which has already ended by today is also a possibility. The model (denoted AS in the following) has the 
following potential : 

F($) = A/4e-^*(Po + (*-$c)'). (62) 

The potential ea. (|52|l has a small minimum in order for the field to be trapped at <I>± = $c + (1 ± \/l — X^Po)/X. If 
X^Pq < 1 the minimum exists and the field plays for a while the role of a quasi-cosmological constant term ($ slows 
down but doesn't oscillate), however if the field has enough kinetic energy it can roll over the barrier. If X^Pq > 1 
there is no minimum and the potential has to be flattened sufficiently for acceleration to occur and fig can reach 
0.55. In section V.B an accurate interval will be given for X^Pq. The only flne tuning in this model is the value of 
which expresses the cosmological coincidence problem : $c is roughly the minimum of the potential and so defines 
the moment when acceleration begins. Once $c is fixed, the beginning of the accelerated stage is given. Apart from 
that, all parameters (M,A,Po,<I>c) take natural values. 

To summarize, we have to take A > 9 in order to satisfy the nucleosynthesis bound, A^Po ~ 1 in order for the 
acceleration to be transient and possibly ending before the present time, and finally <I>c has to be fine tuned so that 
quintessence dominates today. 
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C. Pseudo exponential potential 

For the AS potential ea. (|62|l the minimum, when it exists, is fixed. Hence the acceleration will be eternal for most 
parameter values, see Figure |21 With a straightforward generalisation, allowing a simple coupling between two scalar 
fields, we can obtain a transient acceleration, possibly ending before today, whenever acceleration takes place. The 
introduction of an auxiliary field ^' will control the presence or not of a minimum for The idea is to have initially 
on one hand a minimum of the potential in the $ direction responsible for the acceleration of our universe, and on 
the other hand an evolution of the ^> field such that this minimum disappears in the course of time allowing the 
resumption of matter domination. We will study a two-fields potential of the form : 

V{<f, = M^e-^^(Po + /(*)($ - + (63) 

The AS model [38,] is recovered when / = 1 and g = 0. For the potential ea. (|63|l the minimum is now located at 

The function g {g > 0) can describe a mass term for 4* of the form g oc ^f^, but it is not essential for the dynamics of 
the model. Thus, we will take 5 = for simplicity. 
The minimum Ht)4|) will disappear provided we have 

/(*) < X^Po = /(*c). (65) 

We note in passing that the potential can be rewritten in the form V{^, '9) = Af''/A^e"^*(/(\E'c) + /(*)(A$ - A<i>c)^). 
We will use for / a positive, continuous function which is monotonic in the region 5* > and/or ^E* < 0. As can be seen 
from the condition (|65l) . if / is decreasing, rcsp. increasing, then for smaller, resp. larger, than 'I'c acceleration is 
possible because the minimum ()64|l exists. 

During the evolution of the universe, the field <i> rolls down its potential, which is dominated by the exponential 
part, such that Mpm| ~ M'^m% ^ y ~ ~ M^H-^ with fig = wq = 1/3 during radiation domination, while 

fig = 3/A^, Wq ~ during matter domination until $ approaches its minimum. As long as the field $ is trapped 
at its minimum (|64|l and oscillates around it, the universe undergoes an accelerated expansion with V ^ and 
therefore wq ~ —1 until "if satisfies '5 < resp. ^ > ^'c, if / is increasing, resp. decreasing. At that moment the 
minimum (|64(l disappears, allowing $ to continue to roll freely towards larger values and hence the matter dominated 
regime is resumed. For a given set of parameters (A,Po, ...), the further the initial condition Vl/^ from the longer 
the evolution of '5 toward and so the longer the accelerated regime of our universe. When ^ is initially larger 
(smaller) than '^c, provided / is increasing (decreasing) ^ passes through because dq,V = Af'*e~'^*($ — ^c)^'^ is 
positive (negative), acceleration occurs which is always transient; if on the contrary '^i < "ifc {'i'i > 5'c) quintessence 
domination is not possible. The critical value '^c controls the presence or not of the minimum for <&. Contrary to 
|l6l IT^ where acceleration is always permanent, or |l8j| where acceleration can be either transient or permanent, 
acceleration here is always transient as in . Of course the absence of acceleration is also possible in principle 

for all these potentials but this possibility has to be rejected by observations. 

For the numerical computations presented in section V.C we will use a very simple function / without any additional 
parameter, namely 

/(*) - , (66) 

for which the minimum of <i> diseapears if — ^E'c < ^' < with = X^/Pq. Analogously to the AS model, we have 
to take A > 9 in order to satisfy the nucleosynthesis bound and, once is fixed, <i>c has to be fine tuned in order for 
quintessence domination to occur today. 



D. Pure exponential potential with a varying coupling constant 

Starting from a pure exponential potential V{^) cx e^'^*, we now allow A to depend on the auxiliary field 'i' and 
to vary in time, i.e. we make the generalization A Ae//(^'), and we consider the potential 

y($,*) = Af^e-^''"^*)* . (67) 

A similar idea was used in fss'l where Ae//($) depends on $ in order to create a minimum and to produce an eternal 
acceleration. 
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We will assume that Ag// has a global minimum equal to A in v]/ = 5'„iin = 0. Whatever the precise form of 
Ae//, d,s,\nV = —^eff < for all ^' which imphes that $ is always growing and dq,\nV — —^-^^ is positive, resp. 
negative, if $ and 4" have opposite signs, resp. the same sign. Hence if initially $i < and > 0, resp. '^i < 0, 
(we could take $i > as well with Ae// ~^eff) grows towards positive values while ^' decreases, resp. grows, 
until it reaches, and oscillates around, zero. While <i> < and 5" ~ 0, Ae// A and V e~^*. Hence, as for the 
pure exponential A < 2.3 implies sufficient quintessence domination {^Iq^ ^ 0.55) after the matter dominated stage. 
Once $ > 0, ^' moves away from zero and the function Ag// is again growing, allowing for the resumption of matter 
domination {Qq <^ 1). The current quintessence domination begins when 5* goes to zero and consequently when Ae// 
reaches its minimum. Thus will only depend on A and so will the dynamics of $. The lower A, the slower the 
evolution of $ toward zero and hence the longer the acceleration and/or the quintessence domination regime of our 
universe. Clearly the more remote $i from 0, the stronger this effect. 

The quantity M is solely determined by the conditions to have a realistic evolution of our universe and turns out 
to be roughly equal to the energy scale today M ~ ^ MpHQ. In this potential |Ae//$| < — lnM^/37J^ because 
Xy < 1. So the potential H67(l allows us to use the scaling property of the exponential potential as early as the end of 
inflation without necessarily having kination followed by a stage where the fields are frozen. If Xy^ ^ 1, then after a 
kination regime {wq ~ 1), we have pQ <^ pb and the evolution of the fields is frozen since m| ~ M~'^V ^ 

until V ^ MpH^. For potential H67|l . as for the potential (|63|l . acceleration if it takes place is necessarily transient. 
Moreover the full range < A < 2.3 gives viable models, either with transient acceleration or without acceleration 
because for this potential the domination of quintessence is transient. 

We will take for numerical simulations in section V.D : 

Ae//(f)=A(l+a*2) , (68) 

with a > 0. To summarize, for the potential (|67|l . 5* controls the beginning of quintessence domination while $ 
controls its end. Observational constraints require the following two conditions : A < 2.3 and ~ MpH^. 

V. NUMERICAL RESULTS 

We summarize briefly the observational constraints and the initial conditions used in our numerical calculations. 
We will use the following conservative constraints H41I42I46|I on the quintessence density today r2Q,o, the effective eos 
parameter lUe// and the age of the universe t^: 

0.55 < 17q^o < 0.80, We// < -0.70, iJo^o > 1-330 /i . (69) 

From the equipartition of energy, natural initial conditions suggest that fig.i ^ 10~^ — 10""* When Og^^ < 1 %, 
the parameter window allowed by observations is unchanged and so we will take 

X^^, = = Xy., = 10^2 at = -67 , (70) 

implying ^qa = 3 x 10"'' except for the pure exponential potential. Section V.A, where we take X| - — ■ — 10"'* 
and Xy^ has to be adjusted. 

A. Pure exponential potential 

The exponential potential ea. (|49|) is a viable candidate for quintessence provided the attractor solution is not yet 
reached as explained in Section IV. A, and in particular it can produce transient acceleration. In addition to the 
constraint A > a/2 in order to have transient acceleration (see section IV. A), imposing the observational constraints 
(|69|l only A < 1.975 is allowed. In only the case A < V3 was considered because they used the constraints 
WQfl < —0.60 and ftg.o > 0.60 which lead to a stronger constraint on A. 

Thus in order to have transient acceleration, A is constrained as < A < 1.975 resulting in We// ^ —0.86. Here 
3.7 < X 10"3 < 9.2 with $i = at TVi = -67. 

Different possibilities arise : if (1.82—1.837) < A < 1.837, the acceleration ends before today with a low quintessence 
density ^ 0.62 and Weff ^ —0.75. Surprisingly there is no acceleration at all in the interval 1.838 < A < 1.975 
while 1^(5,0 i$ 0.61 and Weff ^ —0.75. If the observations would constrain Weff so that Weff ^ —0.86, only permanent 
acceleration would be possible (i.e. A < v^). 

In Fig.^ the ratio tend/to is shown as a function of A. 
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FIG. 1: The ratio is displayed as a function of A for the pure exponential potential eg. 14911 . For t > te„d, decelerated 
expansion is recovered {q > 0). All models displayed have a transient acceleration, for those below the dashed line the 
acceleration is ended by today. Initially, at Ni — —67, we take Xy ^ = Vi/pc,i = M^/pc,i = 5.8 x 10~^^^. For this initial value 
Xv,i, the range 1.833 < A < 1.837 yields an accelerated stage which ends before today while the range 1.838 < A < 1.975 
produces no acceleration at all, the range A > 1.975 being excluded by observations, ea. l|69^ . Note that tend +oo when 
\^ V2. 

We should stress that the constraints we have adopted are conservative and a more refined a naly sis would restrain 
the viability of our models. Let us consider the following constraints at la from WMAP alone |26|: 



13.2 Gpc < 3 



dz 



< 14.2 Gpc , 



(71) 



0.12 < ^l„j o/i^ < 0.16 



(72) 



0.67 <h< 0.77 



(73) 



13.1 <to< 13.7 Gyrs 



(74) 



We have finally from the SNIa data at Icr |4£ 



0.26 < n„„a < 0.34 



(75) 



It is easily checked that the range (fSUIl corresponds to ifT^ with the lower, resp. upper bound, corresponding to 
h « 0.81, resp. h « 0.58. However in a more refined analysis the quantities h and flm,o are no longer independent. 



In Fig. El the quantity Dl{z) is plotted for an ever decelerating universe with A = 1.84, Xy^ = 5.8 x 10 
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This model, representative of the scenario without acceleration at all, is characterized by the following cosmological 
parameters ilg.o = 0.583, wq^ = —0.568, Weff = —0.722, Hoto — 0.843. With a correction factor 1.11 < < 1.18, 
the luminosity distances agree with the SNIa data, with uncertainties taken at Icr. Taking the fiducial value — 0.72, 
we would get 0.61 < h < 0.65. For this model we have further J^"^"" — 2.668 so that the constraint (|7T|l yields 
0.56 < h < 0.61 and 0.133 < ^Im.oh^ ^ 0.153 and an age 13.51 < to < 14.72 Gyrs. So this extreme scenario is already 
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FIG. 2: The allowed window (the area between the solid curves) in the parameter space (A, A^Po) is shown for the AS potential 
ea. ib'Zl when the acceleration is transient and all observational constraints, eg. 16911 . are met. The nucleosynthesis bound implies 
A > 9 while transient acceleration implies A < 25. The observations constrain the remaining parameter of the model $c to be 
in the range 237.5 < A4>c ^ 239 for the initial condition "l>i = at A'^; = —67. Above the lower dashed line, qo can be negative 
as well as positive, above the upper dashed line qo > always. For A^Pq values below the window, acceleration is eternal as <1> 
settles at its minimum forever. 

in marginal agreement with the WMAP data alone due to the low value of h, a value h ~ 0.61 being slightly below 
the 2(7 error from the WMAP data alone. 

Let us consider now when WMAP is combined with other CMB data and other data probing the power spectrum 
of the perturbations, we get at la : 



13.7 Gpc < 3 h-^ 



dz 



Hz) 

0.126 < nrn.oH < 0.143 , 



< 14.2 Gpc , 



(76) 
(77) 



0.68 <h< 0.75 , 



(78) 



13.5 <h< 13.9 Gyrs . 



(79) 



We get now the tighter bounds 0.56 < h < 0.58, 0.133 < ilm,oH < 0.142. The value of h is too low and the model 
without any acceleration is clearly in trouble. 

We should however make the following important remark: the uncertainties are obtained from the data assuming 
a ACDM model and a specific model for the perturbations, a constant spectral index for (|71(I - H74(I and a running 
spectral index for (|76|I - H79|) . Our models are not, by definition, ACDM models so that, strictly speaking, the data 
should be processed specifically for each of our models. 



Albrecht & Skordis potential 



For the AS potential ea. (|62|l . only a tiny interval for X^Pq namely 0.985 < A^Pq ^ 1.127 allows a transient 
acceleration if we impose the constraints (|69|l . see Fig. |2 If X^Pq < 0.985, the field <f> stays at its minimum and 
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FIG. 3: The evolution of the densities Qr (dotted line), ^Im (dashed line) and ^Iq (solid line) are shown for the potential eg. lSOII 
with A = 10, Po = 0.164, <l>c = 23.8 ($i = and *i = 5). The quintessence density today is Qq^ ^ 0.661. The age of the 
universe is Hoto — 0.912 and accelerated expansion stops at tend — 0.996 to- The scaling behavior of the quintessence field $, 
ea. lSOL is obvious when the background is radiation or matter dominated, in the past as well as in the future. 





^Q,0 




Hoto 


tn 


WQfi 


90 


1.07 


0.556 


-0.782 


0.838 


0.879 


-0.250 


0.291 


1.06 


0.573 


-0.804 


0.849 


0.907 


-0.303 


0.240 


1.05 


0.592 


-0.828 


0.862 


0.936 


-0.364 


0.176 


1.04 


0.613 


-0.854 


0.878 


0.966 


-0.437 


0.098 


1.031 


0.634 


-0.879 


0.895 


0.997 


-0.516 


0.009 


1.03 


0.637 


-0.882 


0.897 


1.001 


-0.526 


-0.002 


1.02 


0.663 


-0.913 


0.920 


1.042 


-0.634 


-0.131 


1.01 


0.691 


-0.946 


0.948 


1.098 


-0.768 


-0.296 


1. 


0.719 


-0.977 


0.978 


1.194 


-0.914 


-0.486 


0.99 


0.733 


-0.993 


0.995 


1.478 


-0.996 


-0.595 


0.986 


0.732 


-0.994 


0.994 


2.306 


-1. 


-0.598 



TABLE I: Models for the AS potential 1621 with fixed parameters A = 10 and A^c ~ 238.5 are tabulated. All models have 
transient acceleration, for those in the upper part of the Table the expansion is already decelerated today {tend < to). Note 
that the model at the top of the lower part of the Table, A^Po — 1.03, satisfies the necessary condition for decelerated expansion 
today 14311 ■ however this condition is not sufficient as Qq^ — 0.637. 



acceleration last forever while if X^Pq > 1.127 acceleration can occur but flg.o will never reach the value 0.55 because 
the minimum does not exist and the evolution of the field is too fast. If the transient acceleration ends before today, 
again with the constraints (|69|1 the parameter space (A, X^Pq) is even more restricted, see Figure|21 The parameter A is 
constrained to be A > 9 because of the nucleosynthesis bound if the scaling behavior starts earlier than nucleosynthesis, 
while transient acceleration implies A < 25. For this interval, we have 237.5 < A$c ^ 239 when we start with $i = 
at N, = -67. Moreover M ~ IQ-^p^// if we impose = 10""*. 

In the table some examples of a transient acceleration are given with fixed parameters A = 10, A<f>c = 238.5 
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FIG. 4: The evolution of wq (solid line) and the deceleration parameter q (dashed line) are shown for the model of Figure |3 
We obtain here wq,o — —0.491, qo = 0.013 and Weff — —0.874. Again, we note the scaling behavior of wq, eg. 1)50^. In the 
future we have a kination regime {wq ~ 1) and a freezing regime {wq ^ —1) of the field $. 



Po 


^Q,0 


Weff 


Hoto 


tend 

tn 


•WQ,0 


90 


0.17 


0.575 


-0.771 


0.846 


0.889 


-0.207 


0.321 


0.168 


0.600 


-0.801 


0.863 


0.915 


-0.280 


0.248 


0.166 


0.628 


-0.835 


0.885 


0.954 


-0.372 


0.149 


0.164 


0.661 


-0.874 


0.912 


0.996 


-0.491 


0.013 


0.163 


0.679 


-0.894 


0.929 


1.02 


-0.564 


-0.075 


0.162 


0.699 


-0.916 


0.949 


1.05 


-0.649 


-0.181 


0.16 


0.740 


-0.960 


0.995 


1.15 


-0.850 


-0.444 


0.15 


0.746 


-0.989 


1.008 


2.49 


-0.993 


-0.611 


0.14 


0.716 


-0.990 


0.978 


3.62 


-0.995 


-0.569 


0.13 


0.681 


-0.990 


0.947 


5.02 


-0.996 


-0.518 


0.12 


0.643 


-0.989 


0.917 


6.80 


-0.998 


-0.463 


0.11 


0.602 


-0.988 


0.889 


9.06 


-0.998 


-0.401 


0.1 


0.559 


-0.987 


0.862 


11.98 


-0.997 


-0.336 



TABLE II: Models for the potential l|8()|l with fixed parameters A = 10 and A$c ~ 238 are tabulated. The four upper models of 
the table have a transient acceleration already ended by today, see Fig. |S| The model at the top of the lower part of the Table, 
Po = 0.163, satisfies the necessary condition for decelerated expansion today H43|l . however this condition is not sufficient as 
^Q,o = 0.679. 

and we vary Pq. All the points for which the transient acceleration is ended by today have Weff ^ —0.88. If Weff 
is constrained from observations to satisfy Weff < —0.88, acceleration ended by today is not possible in the AS 
potential. Nevertheless transient acceleration ended by today with the AS potential is easier to achieve than any 
transient acceleration using the pure exponential with We// ^ —0.86. 

Whatever A and Pq, the maximum value for tend /to is 3 because Pq is bounded from below (see also Figure [2Jl. 

At the end of Section V.C, a model representative of all four models with acceleration ended by today will be 
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FIG. 5: The ratio where fe„d is the time at which accelerated expansion ends, is shown for the potential eg. (1801 versus the 

ratio with ~ \^/Po. The displayed models have the fixed parameters A — 10, "I>c = 23.8 with the initial conditions 

— and 'i'i = 5, while we vary the parameter "ifc, or equivalently, Pq. All models displayed have transient acceleration but 
only those below the dashed line have a decelerated expansion today. The range of values 'I'o/'I'c < 0.95 is excluded by the 
observations for this set of parameter values. 

analysed in the light of the constraints H71|l - H79() . 

C. Pseudo exponential potential 

We can illustrate the considerations concerning the potential ea. Hfci3|) (with g = 0) using the function defined in 
ea. (|(j6|l . The potential reads then 

V^—e-^'^{^'l + -i'^{X<P-X<Pcf), (80) 

and the minimum for $ disappears if — ^l/c < Vl/ < ^'c with "ifc = X^/Pq. Here we will use the following initial conditions 
$i = and 'i'i = 5. Once <i>i and '^i are given, we can choose ^I^c so that < ^I'c < ^'i while it is always possible 
to find an appropriate <i>c. In a way analogous to the numerical calculations with the AS potential, we keep the 
parameters A = 10, A'l'c = 238 fixed and we vary '^c- Also, as for the AS model, we have M ~ 10^^/9^/ if we impose 
X^,, = 10-^ 

In Fig. O the evolution of densities are plotted for the parameter Pq = 0.164 (5*^ ~ 4.05). This set of parameters 
induces a transient acceleration which ends before today. In the scaling regime, Qq ~ 0.04 during radiation domination 
and Qq ~ 0.03 during matter domination, while ^Iq.q ~ 0.661. 

In Fig.0I the eos parameter is plotted with wq^ ~ —0.491 and We// = —0.874 and the scaling behavior is evident: 
WQ ~ 1/3 during radiation domination, wq ~ during matter domination. When the field $ settles at its minimum, 
WQ ~ — 1 and this stage is followed by a stage of kination and freezing out of $ once the minimum is left. The 
deceleration parameter q is also plotted implying an actual deceleration as go — 0.013. Acceleration {q < 0) begins 
at z ~ 0.658 and finishes at 2 ~ 0.0035 when the age of the universe is iend/^o — 0.996 while the present age of the 
universe satisfies H(^tQ ~ 0.912. 

In Fig. El the ratio is plotted as a function of the ratio ^'o/^'c- Clearly, the more remote {'i'o) from the 
longer the acceleration regime. As soon as the minimum disappears, $ stops oscillating and continues to roll down 
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FIG. 6: The evolution of the quantities wq (solid line) and q (dashed line) are shown versus redshift for the po- 
tential ea. l|8()|l with the fixed parameters A = 10, $c = 23.8 while Po takes the values, from bottom to top, 
0.160,0.162,0.163,0.164,0.166,0.168,0.170 (<l>i = and *i = 5), see Table |ll| Note that the third model from bottom, 
Po = 0.163, satisfies the necessary condition for decelerated expansion today (143 l l. however this is not sufficient because 
^Qfi > 0.55. The models with Po = 0.164, 0.166, 0.168, 0.170 produce a transient acceleration which ends before today (go > 0). 

its potential. Contrary to the AS potential tend /to can be as large as we want because Po, or equivalently ^fc, is not 
bounded from below. 

In Fig. El Wq and q are plotted against redshift: the cases Pq = 0.164,0.166,0.168,0.170 induce an acceleration 
ended by today. Viable models where the acceleration ends before today have Weff ^ —0.88. 

In the table |n] some examples for this model are shown with $i = 0, vE'i = 5 and again A = 10, <I>c = 23.8 and we 
vary Pq (or equivalently ^'c). 

In Fig.O the quantity is plotted for a model for which the acceleration is already finished (A — 10, Pq = 0.164, 
$c = 23.8, $i = 0, = 5, Qm,o = 0.339, Hoto = 0.912). With a correction factor 1.03 < ga < 1.10, the luminosity 
distances agree with the SNIa data, with uncertainties taken at la. Taking the fiducial value h\ = 0.72, we would 
get 0.65 <h< 0.70. 

We can repeat for this model a more refined analysis of the kind done for the other model (pure exponential) shown 
in Figures m Eni see Section V.A. We have J^""" = 2.992 so that the constraint ^ yields 0.63 < h < 0.68, 
0.135 < flmfih^ < 0.157 and an age 13.11 < to < 14.11 Gyrs. So this model is in good agreement with the 
WMAP data alone and the SNIa data. If we consider the tighter constraints (|76|l . we obtain 0.63 < h < 0.66, 
0.135 < rim.oh'^ < 0.145 and an age 13.51 < to < 14.11 Gyrs, stiU allowed by (f7H|l - (f77jl . 

D. Pure exponential potential with a varying coupling constant 

We now consider the potential ea. (|67l68|l 

V = exp(-A(l + 0*2)$) , (81) 

for which quintessence domination occurs when A < 2.3 and acceleration today takes place provided M ~ ^jMpH^- 
We have taken a\ = 12 AA with = —5 and '^i = 2 and initially Xy^ = 10"**. Thus the dynamics of the fields starts 
initially, that is, there is no freezing out {V ~ H^). 
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FIG. 7: The evolution of the quintessence density is shown for the potential eg. 18111 and same models as in Table lITTIl with, 
from top to bottom, A — 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.85. For the models A = 1.8, 1.85 (A > %/3), there is an intermediate 
stage where Qq J7lb (cfr ea. lSOH l starting around A'' ~ —1.4 {z ~ 3); for the models A = 1.2, 1.3, 1.4, 1.5, 1.6, 1.7 (A < \/3), 
there is an intermediate stage where Qq 1 (cfr ea. l51ll l. For all models displayed here, Qq —> asymptotically. Note that 
for a — (pure exponential), the models with A — 1.2, 1.3, 1.4 (A < y/2, see ea. 1151^ '1. would yield eternal acceleration and 
domination of the quintessence energy density. We see that the smaller A, the longer quintessence dominates. 
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-0.715 
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0.755 


-0.792 


0.973 


56.31 


-0.691 


-0.282 


1.4 


0.717 


-0.778 


0.937 


63.94 


-0.669 


-0.220 


1.5 


0.680 


-0.765 


0.907 


3.62 


-0.649 


-0.161 


1.6 
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TABLE III: Models for the potential ea. (l8H with fixed parameter values aX = 12.44, and initial conditions <l>i — —5, — 2, are 
shown here. All models have transient acceleration except for the model at the bottom of the Table which has no acceleration 
at all. Note that the necessary condition 1431 for decelerated expansion today is satisfied for A = 1.8 but it is not sufficient as 
nQ,o = 0.580. 

In Fig. [71 the quintessence density evolution is plotted for different values A = 1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.85. 
The effect of the value of A is shown where, for fixed initial conditions, the models today are roughly equivalent, but 
the duration of the transient acceleration depends strongly on the value of A. 

In Fig.lSI wq and q are shown and we note a permanently decelerating universe for A — 1.85, however all the models 
with no acceleration at all have Weff ^ —0.74 and are marginally viable. Note that the oscillations of the auxiliary 
field ^ around z ^ 1 — 3 translate into oscillations in wq, but it is impossible to see them in the luminosity distance 

dL 

In the table mil below some examples for this model are shown with <I>i — —5, '^i — 2, aX = 12.44 and we vary A. 
Scenarios with no acceleration at all in model D are of the same type as for the pure exponential, model A. The 
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FIG. 8: The evolution of wq (solid line) and of the deceleration parameter q (dashed line) is shown for the potential eg. lSlH 
with the same parameters as in Figure 7, and with growing values of A from bottom to top. For the model A = 1.85, there is 
no accelerated expansion at all. Note that wq exhibits oscillations in the range 2: ~ 1 — 3 because the auxiliary field 'J reaches, 
and oscillates around, 0. 



analysis performed at the end of Section V.A is therefore representative for model D. 



Models 


Scenarios 


Parameters 


Fine tuning 


Pure exponential 
V = M*e-^* 


eternal acceleration 


< A < \/2 




transient acceleration 


^< A< 1.837 


no acceleration at all 


1.838 < A < 1.975 


AS 

V = M^e-^* [Po + (* - $c)^] 


eternal acceleration 


9 < A 

A^Po < 0.985 


<l>c depends on i.e. 


transient acceleration 


9 < A < 25 
0.985 < A^Po < 1.127 


AS generalisation 

V = M*e-^* [Po + - ^.f\ 


transient acceleration 


9 < A 

A^Po < *o < 


<l>c depends on i.e. 


Exponential with a varying A 

y ^ jy^4g-A(l+Q*2)<l. 


transient acceleration 


< A < 1.84 




no acceleration at all 


1.84 < A < 2 



TABLE IV: The possible scenarios for our one-field models (upper part of the Table) and our Double Quintessence models 
(lower part) are summarized. As can be seen, for the one-field models, eternal acceleration covers most of the allowed parameter 
space. Note that for the Double Quintessence models, the acceleration is necessarily transient. All four models can produce 
scenarios satisfying 16911 with a transient acceleration already ended at the present time. 



VI. CONCLUSION 

Recent observational data suggest that some unknown component, called Dark Energy, contribute to about two 
thirds of the present total energy density filling our universe. Though the accuracy of the existing observations allow 
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FIG. 9: The following models are plotted: the SNIa data with la errors corresponding to flat A-models with SIa.o = 0.74 (upper 
solid line), i}A,o = 0.66 (lower solid line), the Einstein-de Sitter universe (dotted line), the potential eg. 1)80^ for Po — 0.164 of 
Table HTl and accelerated expansion ended by today (long-dashed line), the pure exponential potential ea. l|49^ for parameter 
values A — 1.84, — 5.8 x 10"^^'^ and no acceleration at all (short-dashed line). The luminosity distances dhiz) for a given 
model with Hq can be compared by varying ba, e/^DL{z) = ^Dl{z) — Ho,AdL{z), where //o,a (and the corresponding /ia) 
is some fixed fiducial value. The correction factor 1.03 < ca < 1.10 will bring the first model (pseudo exponential) inside the 
SNIa data {la uncertainties), 1.11 < ca < 1-18 is needed for the pure exponential. 

to constrain already at the present stage to some extent the possible Dark Energy candidates, a very large number of 
models are still permitted. It is therefore interesting to investigate all possible scenarios and we have also investigated 
in this work specific models where the Dark Energy sector is made of two coupled scalar fields. We have studied 
numerically two one-field models, the pure exponential in Section V.A and the AS model in Section V.B, and two 
Double Quintessence models, which can be seen as extensions of the corresponding one-field model, in Section V.C and 
V.D. We were interested in particular in scenarios for which the (recent) stage of accelerated expansion is transient. 

Constraining the models with the observations, we have found the allowed window in the corresponding parameter 
space of each model. Investigation of these models has revealed that the following three possibilities can arise: the 
present acceleration is transient and still going on; some accelerated expansion did take place in the recent past but is 
already finished by today; finally no acceleration at all, this latter possibility being marginal. The two first scenarios, 
i.e. a transient acceleration either finished or not by today, can be obtained in all four models. A scenario with no 
acceleration at all is obtained only for the two one-field models V.A and V.D provided the cosmological parameters 
take their values at the edge of the allowed range, ^lm,o ^ 0.4 and h < 0.65. We summarize the possible dynamics of 
our models in Table Hvl 

In both Double Quintessence models studied in Sections V.C and V.D, for the allowed window in parameter 
space where accelerated expansion takes place, the dynamics of the model is such that the acceleration is necessarily 
transient, hence for these two models acceleration if it takes place will eventually come to an end either before today or 
in the future. In the model of Section V.C, in addition to the tracking field an auxiliary field ^ is introduced which 
controls the presence of a minimum for the $ field and the duration of the Universe transient acceleration. In the model 
of Section V.D, it is the auxiliary field 4", and hence Xeff, which induces accelerated expansion when it reaches its 
minimum. In contrast to our Double Quintessence models where acceleration is necessarily transient, for the one-field 
models we have studied the acceleration is typically eternal and covers most of the allowed parameter space. In view 
of the theoretical problems posed by eternal acceleration, all viable scenarios with transient acceleration constitute a 
wellcome alternative. On the other hand, the cosmological coincidence is not solved here and requires some amount 
of fine tuning on one of the free parameters. We insist that all the scenarios studied here are in agreement with 
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FIG. 10: The relative error for the luminosity distance dL{z) with respect to the A model with Oa.o = 0.66 is shown. One has 
- 1. The following models are plotted: the Einstein-de Sitter universe with ba = 1 (dotted hne), the 
0.164 of Table im and accelerated expansion ended by today, with ga = 1 (lower long-dashed line) and 



dL(h) 



SA 



<1l, A{hA) 

potential eg. 11801 for Po 

ea = 1.03 (upper long-dashed line), the pure exponential potential eg. 1491 for parameter values A = 1.84 and Xy i = 5.8 x 10"^^'^ 
and no acceleration at all, with ea = 1 (lower short-dashed line) and ca = 1.11 (upper short-dashed line). As seen from the 
Figure, by varying ca, and ca > 1, we improve agreement with the SNIa data, except for small z < 0.3 where the relative 
departure can be large. Clearly, this is impossible with the Einstein-de Sitter model.x 



observations, in particular with the Hubble diagram H{z), or the luminosity distances d^i^z) as a function of redshift, 
as reconstructed from the Supernovae data leading to the possible interpretation of a flat universe with f^A.o — 0.72. 

As mentioned in the Introduction, it is well known that many-fields inflationary models can produce primordial 
fluctiations spectra with a characteristic scale. The question arises naturally whether similar effects can be produced 
here. In fact, the luminosity distances can exhibit a characteristic scale if the equation of state of dark energy 
undergoes a phase-transition. Such cases were considered in [50l |. where the quantity w{z) was taken with a step-like 
structure. We have checked numerically that such models could indeed exhibit a characteristic scale in their luminosity 
distances. In our Double Quintessence models where accelerated expansion is ended at the present time, we have also 
a large variation of w(z) at low redshifts, typically from w{z) ~ —1 up to some higher value. However, as can be 
seen from Figures El and |S1 this variation takes place at very low redshifts < z < 0.3 and, though significant, this 
variation is not sharp enough. This is the reason why no characteristic scale is seen in the corresponding luminosity 
distances. Actually, it was already emphasized that the luminosity distances are not very sensitive to large, smooth, 
variations of the equation of state, or equivalently of the quantity t«(z), at low redshifts z. Results obtained here are 
a particular illustration of this property. 

In the light of our results, it is clear that eternal acceleration could be challenged in two ways, either the acceleration 
is transient and will end at some time in the future, either it has already ended by today. Results obtained here show 
that the second possibility must be taken seriously if the observations allow a rather high matter content today, viz. 
^m,a ^ 0-35 and also large variations of the eos parameter w{z) at low redshifts. In particular, if it turns out that 
the data require wq^ < —0.5, this possibility is ruled out. 
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